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Abstract 

We study the self-trapping of quasiparticles (electrons, holes, exci- 
tons, etc) in a molecular chain with the structure of a ring, taking into 
account the electron-phonon interaction and the radial and tangential 
deformations of the chain. A discrete system of equations is obtained 
and solved numerically. The analytical solutions for the wave function 
of a quasiparticlc and for the molecule displacements that determine 
the distortion of the ring, are also obtained and solved in the contin- 
uum approximation. The numerical solutions of the system of discrete 
nonlinear equations reveals several regimes of quasiparticle localisation 
in the chain which depend on the values of the parameters of the sys- 
tem. It is shown that the transversal deformation of the chain favours 
the formation of a soliton. 



1 Introduction. 

The study of various physical properties of carbon nanotubes has received 
a lot of attention during the last few years. In particular, B. Hartmann et 
al [1] have recently studied the electron-phonon interaction on a hexagonal 
lattice. Their approach was, however, limited in its scope as it did not take 
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into account the possible displacement of atoms in the direction perpendic- 
ular to the nanotube. To address this problem, we have to introduce three 
phonon fields, two of which are tangent to the tube while the third one is 
perpendicular to it. The introduction of such a transversal field is not very 
straightforward. It requires the generalization of the phonon Hamiltonian to 
describe the bending elasticity of the tube. Such a generalization is presented 
in this paper for the case of a nanocircle. A study of the electron-phonon 
interaction in a simple model where such a transversal field is introduced, 
paves the way towards the study of a more relevant model describing car- 
bon nanotubes. It also describes various physical properties, including the 
conductivity and conformational states, of other physical systems such as 
circular macromolecules, cyclic organic, inorganic and semi-inorganic poly- 
meric macromolecules. 

Our model describes a circular closed chain of monomers which are al- 
lowed to move in the plane spanned by the circle at rest. Thus, it is a 
two-dimensional model where deformations in the direction perpendicular 
to the circle plane, are not allowed. One can think of this model as a crude 
cross section of a carbon nanotube, though obviously, the hexagonal symme- 
try of the carbon nanotube, and its three-dimensional nature require a more 
general description. Nevertheless, our hope is to gain some understanding 
of how one can take into account the transversal displacement of atoms. 

At the same time, our model does describe real physical low-dimensional 
cyclic nanosystems that include a large class of end-linked natural and 
synthetic low-dimensional macromolecules and nanosystems with a circu- 
lar structure. Amongst them one finds the cyclic biological macromolecules, 
such as intact DNA from most prokaryotes and viruses, mitochondria plas- 
mid DNA [2, 3, 4], cyclic DNA macroarrays used for the systematic screening 
of gene expression [5, 6, 7]. Large circular single stranded RNAs form ge- 
nomic material of viroids [8, 9]. Some of the oligoribonucleotides with known 
biological activity are also cyclic molecules [10, 11, 12, 13]. The so called 
"second messengers" also have a circular structure. Amongst them are the 
messengers which regulate the channels involved in the initiation of impulses 
in neurons responding to odors and light, messengers which function as 
calcium regulated transcription factors or as a substrate for depolarization- 
activated calcium calmodulin-dependent protein kinases I and II (e.g., cAMP 
Response Element-Binding Protein, known also as CREB protein) [14]. In 
addition, there are also messengers which induce and control transcription 
from the proximal promoter of the human aromatase responsible for the 
metabolism of steroids [15]. Therefore, our study of distortions of molecular 
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nano-circles due to electron-phonon coupling could bring some insight into 
the structure- function relationship of biological macromolecules. 

End-linked organic, inorganic and semi-inorganic polymers are now known 
as functional materials with applications as biologically active species, cat- 
alyst supports, charge carriers etc. For instance, the long-chain circular 
esters are used for soft lenses, cyclic trimeric phosphazene rings are valuable 
technological materials for their high heat resistance and flame retardant 
behaviour. Polysiloxanes, which are the main components of silicone prod- 
ucts, are one of the most important cyclic materials ([16] and references 
cited therein). 

The circular structure of these macromolecules, their relative softness, 
deformability and ability to form cross-linking suggest modelling such sys- 
tems by a nanocircle of periodically placed monomers of a given length. 
Under certain conditions one or several monomers of the nanocircle can be 
excited by an external field (light or electromagnetic radiation, trapping of 
extra electron or hole, etc.). Such an excitation, called in what follows, a 
quasiparticle, can move onto the neighbouring monomers and cause a local 
distortion of the ring due to the electron-phonon coupling. In a linear chain 
it can lead to the change of properties of the quasiparticle, depending on the 
strength of the coupling. At weak coupling the quasiparticle remains almost 
free, at strong coupling it is trapped on one site (a small polaron is formed), 
while at intermediate values of the coupling a self-trapping occurs within 
few lattice sites and a nonlinear soliton-like state is formed [17, 18, 19]. It 
turns out that the physical properties of systems at different regimes of lo- 
calisation differ significantly, and, therefore, the classification of systems by 
their electron ground states is important. 

One can expect that, depending on the strength of the electron-phonon 
coupling, similar regimes take place in a very large end-linked circle. How- 
ever, this is not clear for a circle of a finite and not very large length. This 
problem is studied in the present paper. We start by deriving a Hamiltonian 
describing the elastic properties of a nanocircle and a Hamiltonian describing 
the interaction between the elecron and the phonon fields. We then derive 
the corresponding system of discrete equations and compute solutions for 
this model. We show that the model exhibits several classes of solutions 
depending on the values of the parameters. In particular, we show that the 
electron-phonon interaction leads to the self-trapping of a quasiparticle and 
introduces a deformation of the chain. 
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2 The nano-circle 



We consider a molecular chain which consists of periodically repeated units 
(atoms or groups of atoms) bound by chemical or hydrogen bonds. We as- 
sume that these bonds between the units are directed in such a way that, 
for a molecular chain with free ends, the minimum of the free energy corre- 
sponds to a linear one-dimensional structure in which chain sites are placed 
at the equilibrium distance ao and all intersite bonds are directed along the 
chain axis. 

We focus our attention on a study of the molecular chain which is closed 
to form a circle (before any deformations are allowed). The equilibrium 
positions of the chain sites are determined by the radius- vectors 

R^ = Ro[e x sin(na + Oq) + e y cos(na + #o)] (1) 

where the centre of the coordinate system is placed at the centre of the circle, 
Rq is the circle radius, n = 1, . . . , N enumerates the chain sites, a = 2ir/N 
where N is the number of sites in the ring, and Oq is an arbitrary angle. The 
equilibrium distance between the sites in the molecular circle is 

D { Z +1 = = " 4 0) l = 2i?osin(|) = a (2) 

which, in general, can differ from ao- Thus, the radius of the circle is Rq = 
a/(2sin(^)). Note that if there is a bending of chemical bonds at each 

site n, this bending can be described by the angle (3 n °^ between the vectors 
n(°) - _ an A n(°) _ _ 

C ° S = inV^llnro'r 1 ! = C ° S(a) ' (3) 

i.e., /?i 0) = a. 

Due to the displacements the real positions of the sites are 

R n = R^+r n , (4) 

where r n is the displacement of the site n from the equilibrium position. 
Note that we can represent the vector r n using the following orthogonal 
unit vectors eh/ 1 and en 

effl = g x cos(na + Oq) — e y sin(na + 9q), (5) 
effl = e x sin(na + Oq) + e y cos(na + Oq). (6) 
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We then have 

where r tj „ = u n is the displacement tangential to the circle while r r>n = s n is 
the displacement perpendicular to the circle (the tangential and radial com- 
ponents of displacement respectively). Thus the real distances between the 
sites and the bending angles differ from their equilibrium values: Z) n>n +i = 
a + d n ,n+i and [5 n = a + A/3„. 

The vector, which links the site n with the n + 1 one, is _D„ iri +i = 

R n+ i - R n = £>n} l+1 + d n ,n+i where d n , n +i = r n +i - r n . So, the distance is 
given by 

n in I i"^ \2~ , (-^n,n+l ' d n ,n+l) 

= |Ai,n+l| = \J{D n ' n+1 +d ntn+ iY W OH ; = a+(7„, 

(8) 

where 

= cos(-)(u„ + i - u n ) + sin(-)(s„ + s n+ i). (9) 

To obtain this expression we have assumed that the relative site dis- 
placements are small, i.e. <C a 2 . In this case the deviations of the 
bending angles from their equilibrium values are also small, and we can write 

cos(/3 n ) = cos(a + A/3 n ) (10) 
= cos(a) cos(A/3 n ) — sin (a) sin(A/? n ) « cos(q) — sin(a)(A/3 ra ). 

For cos(/? n ) we can put 

{Dn— l,n ° -Dji,n+l) 



cos/3„ 



l-^n— l,n| \D nyn J r \ | 



(Dn-l,n - + (^i-l,n ' d n ,n+l) + (A^i+l ' d n -l,n) . . 

(a + C/ n _!) (a + C/ n ) 1 j 

, . a u n+ i - u n _i a 2s n - s n+ \ - s n _i 

cos a — sin a sin h cos 



2 a 2 a 

Thus, we have arrived at the following expression for A/3 n : 

AP n = ^, (12) 

where 

OL OL 

A n = sin -{u n+ i - Un-i) + cos -{2s n - s n+ i - s n _i). (13) 



5 



3 The Hamiltonian 



Next, we assume that the individual isolated monomers are characterised 
by nondegenerate electronic excitations with energy E$. In the approxi- 
mation of the nearest neighbour interaction, the total Hamiltonian which 
describes the electronic excitations (or, in the tight binding approximation, 
the additional electrons) and harmonic vibrations in a molecular chain can 
be written as 

Jn,n+1 (Bt n B n+1 + B* n+1 B n j)+H ph , (14) 

n 

where B^, B n are the creation and annihilation operators of the quasiparti- 
cle at the n-th site of the circle, £ n = Eq + V n is the on-site electron energy, 
where V n is the change of the electron energy due to the influence of the 
neighbouring sites. Furthermore, J n ,n+i describes the quasiparticle trans- 
fer between the nearest neighbours, and H p h is the Hamiltonian of the site 
displacements. The electron-phonon interaction is determined by the depen- 
dence of V n and J n ,n+\ on the positions of the neighbours of the n-th site. 
The term J n ,n+i depends on the distance between the sites D n ^ n+ \ = a + U n 
and, within the linear approximation with respect to the displacements, we 
can put it as 

Jn,n+1 = J — G2U n , (15) 

where J = J (a) and G2 = —dJ /da. 

V n is determined by the influence of the neighbours; thus it depends on 
the distance between them and on the bending angle of the chemical bonds: 
V n = V n -i yn (a + U n -i) + V n , n +i(a + U n ) + V n (/3 n ). Hence the on-site electron 
energy, within the displacements' linear approximation, can be written as 

+ XI (U n + tfn-l) + XlK 

OL QL QL Ot 

= £0 + (Xi cos - + X2 sin -)(u n+1 - u„_i) + 2(xi sin - + \2 cos -)s n 
a a 

+(Xi sin - - X2 cos -)(s n+ i + s n _i), (16) 

where xi = dV n ± l , n {a) / da, X2 = (1/ 'a)dV n (f3 n ) / 'df3 n \p n=a . The parameters 
Xi, X2 and G2 determine the interaction of the quasiparticle with the chain 
deformation (exciton/electron-phonon coupling) in the linear approximation 
relative to the site displacements from their equilibrium positions. 
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Due to the symmetry we see that V n ((3 n ) = V n (—(3 n ), and we can set 
V n (0n) = V n (cos{P n )). Therefore, 

X2 = (l/a) d ^^\p n=a = -(l/a)F>n(«). (17) 

Thus, \2 ~ sin(a) and X2 = at a = and in the linear chain the expansion 
of V n ((3 n ) in series of A/3 n begins with the quadratic term while the linear 
term vanishes. 

The Hamiltonian of the site displacements, H p h, consists of the kinetic 
and potential energies of displacements. The potential energy consists of 
the potential energies of the pairwise interactions, which depend on the 
distance between the sites, and on the bending contribution. Therefore, the 
total Hamiltonian of the phonons takes the form: 

1 ~^ 

H p1=2^(m + K(\Rn+i-Rn\-a ) 2 + K b a 2 ((3 n -(3 ) 2 ). (18) 

Here M is the mass of a lattice site, p n are the momentum operators canon- 
ically conjugate to the operators of displacements f n , k is the elasticity 
coefficient due to the change of distance between the molecules from their 
equilibrium value ao, and K b is the elasticity constant responsible for the 
bending of the bond (for small bending angles (3 n ). 

As our chain is closed to form a circle, we have /3 n = a + A/? n = a + 
aA n . From the definition of A n we see that ^2 n A n = and, therefore, the 
equilibrium distance between the lattice sites on the circle remains the same: 
a = ciq. Therefore, the total Hamiltonian for phonons in a circle takes the 
form: 

n 

+ K b (sm^(u n+ i - u n -i) +cos^(2s„ - s n+1 - s„_i)) 2 ) + W h 

= \T.{% + % + KU 2 n + Kb A 2 n )+m, (19) 

n 

where p n and q n are the momentum operators canonically conjugate to the 
operators of the tangential and radial components of the site displacements, 
respectively, and W b = (l/2)Ki,a 2 a 2 N = 2ir 2 n b a 2 /N is the energy required 
to bend the linear chain into a circle (which we will drop from now on as it 
is a constant). 
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Note that both U n and A n are invariant under the global translations of 
the circle and so is the Hamiltonian (19), as one would expect. 
Thus, the total Hamiltonian is given by 

H = H e + H ph + H int , (20) 

where 

H e = ]T (Sq B n B n - J (BlB n+1 + Bl +1 B n )) , (21) 

n 

H in t = E ( B l B n[Xl (Un + U n -i) + X 2 A n ] + G 2 (B n B n+ l + Bj +1 B n ) U n ) , 

n 

(22) 

and H pn is given by (19). 

These expressions are written in the site representation. We can trans- 
form them into the quasimomentum representation by performing the uni- 
tary transformation 

B - = ^= (23) 



S « = 7^I>^ ^ = E^^, (24) 

where fc( or g) = 27ri//iV, i/ = 0, ±1, ±(iV/2 - 1), N/2 for iV even and v = 
0, ±1, ±(N — l)/2 for iV odd. This allows us to rewrite the quasiparticle 
Hamiltonian (21) and the phonon Hamiltonian (19) as 

H e = J2m B lB k , (25) 
k 

q fj, 

(26) 

where fj, = l,r, and 

£(k) =£ - 2J cos(fc), (27) 
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= 4^ S m 2 (|)[co S 2 (|)+4 CS in 2 (|)co S 2 (|)], 
4i<l) = 4^[ S m 2 (f)co S 2 (|) + 4 C co S 2 (|)sin 4 (|)], 

f(q) = ^sin(a)sin( g )[l+4csin 2 (|)]. (28) 

Here c = k^/k. 

Applying the unitary canonical transformation 

U hq = cos(6(q))Q hq + i sm(6(q))Q 2 , q , U hq = cos{6(q))P hq - i sm(6(q))P 2 , q , 
U ryq = i sin{6(q))Qi A + cos(6(q))Q 2 , q , U r>q = -i sm(6(q))P hq + cos(0(q))P 2 , q , 

(29) 

in which 

tan(2%)) = , A{q) = u f {q) _ J2 {q)i (30) 

we obtain the phonon Hamiltonian in the diagonal form 

H * = ^E(4^+^ 2 w^L^)' ( 31 ) 



2 

3,1 



where j = 1,2 and 



tij(q) = i( W 2 ( g )+^ 2 ( g )-(-l)^(c, 2 (^)-^ 2 ( (7 ))co S (20(g))+2/(g) S in(20(g))]). 

(32) 

Then we can introduce the operators of creation }y ■ and annihilation 
b q j of phonons by performing a further canonical transformation 



^ = V2t^) (6w " + 6L « ) ' ^ 



^ = -<\/^(^--6U,-). (34) 



2 

and find that the total Hamiltonian can be rewritten as 



H = J2S(k)BlB k + J2 m M(i,A 3 + (35) 

k j,q 

+-7^ E q)B\ +q B k {b q , + &!„■) (36) 
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where the coupling functions $j(k,q) are given by 
with 

F 1 (k,q) = Ai[F l (k,q)cos(6(q))+F r (k,q)sm(6(q))}, 

F 2 (k,q) = A[F r (k,q)co S (e(q))-F l (k,q)sm(8(q))}. (38) 

Here 

Fi(k,q) = [( Xl cos(|) + X2sin(|))cos(|) + G 2 cos(|)cos(fc + |)]sin(|), 

/j \ . ,a 2 q a. . 2 .q. . a. .. q. ,q. 

F r (k, q) = Xi sm(-) cos (-) + \2 cos(-) sin (-) + G 2 sm(-) cos(fc + -) cos(-). 



Note that $>*{k,q) = $j(k + q,-q). 

At large values of N the parameter a is small (N — ► oo, a — > 0). At 
a = we have a linear molecular chain. In this case, according to (28), 
f(q) = and the longitudinal Qi :Q = Ui >q (along the chain axis) and the 
transversal Q 2 ,q = U T) q displacements do not mix (9(q) = in (29)). A 
linear chain is characterised by the longitudinal acoustic waves U\ A with 
frequency 

u>?(q) =4^sin 2 (|), (40) 
and the bending waves U r>q with frequency 

= 16^sin 4 (|). (41) 

In a linear chain in the linear approximation with regard to the displace- 
ments, quasiparticles interact only with the longitudinal acoustic phonons 
because 

F 1 (k,q) = m(k,q) = 4i[( X i cos(|) + G 2 cos(k + |)] sin(|) (42) 

and F 2 (k,q) = AF r (kq) = 4x2sin 2 (|) = (as \2 ~ sin(a) and so \2 = at 
a = 0). 

Note that in a long linear chain, N which is included in the definition 
of the one-dimensional quasimomentum q, determines the main region of 
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quantization. So, we have q = 2ttv/N = au, v = 0, ±1,.... On a circle 
the longitudinal (tangential) and the transversal (radial) displacements are 
mixed in a way determined by the angle 9{q). From the definition (30) we 
see that tan(2%)) = - tan(20(-g)). This means that 20(q) = -20(-q) or 
26 (q) = ir — 26(— q). We choose the first relation and, consequently, 6(0) = 
and -tt/2 < 20(q) < vr/2. 

According to (32), in a circle with N sites, the phonon Hamiltonian has 
two normal branches with different q which take N values, i.e. 2N frequen- 
cies. Note that in a circle with iV particles, each of which has two degrees of 
freedom, the number of normal vibrations is 2N — 3 because there are two 
degrees of freedom which correspond to the motion of the centre of mass of 
the system and one degree corresponds to the rotation of the circle. And 
indeed, there are three zero frequencies in the vibrational spectrum (32), 
namely: fii(0) = Q 2 (±2-ir/N) =0. At q = 0, 0(0) = and th e nor mal 
coordinate Qi jq =o corresponds to the site displacements u n = 1/y/MNQi^ 
and s n = 0, i.e. this is the rotation of the circle. At q = ±2tt/N = ±a 
0(±a) = ±7r/4, the normal coordinates Q 2) ± a describe the site displace- 
ments: f n = e x ^2/(MN)Qi m2>a and r n = e y ^2/(MN)Q re2:a and cor- 
respond to the motion of the centre of mass in two mutually orthogonal 
directions. Here we introduce the notations Q re j,a = [Qj,a + Qj, a ]/2 and 
Qimj,a = i[Qj,a~ Qj a ]/% after taking into account the relation Qj- a = Q*j,a- 
It is natural that the quasiparticles do not interact with these normal coor- 
dinates. From (38) and (39) we see that Fi(/c,0) = F 2 (k,±a) = 0. 

For the values q = and q = ±a we have two nonzero frequencies which 
have no analogue in an open chain. The frequency 

nftO) =4^sin 2 (|) (43) 

is the frequency of the total symmetrical vibration which corresponds to the 
circle breathing, i.e. to the deviation of the circle radius from its equilibrium 
value Ro, because the normal coordinate Q 2a =q describes the site displace- 
ments s n = l/\/MNQ 2 fi and u n = 0. The doubly degenerate frequency 

n?(±a) = 8^sin 2 (|)cos 2 (|)[l + 4csin 2 (|)] (44) 

is the frequency of vibrations which correspond to the nonsymmetrical de- 
formation of the circle into two mutually orthogonal directions. The normal 
coordinates Qi,± a describe the site displacements fn^ = &n) Aj (j = 1,2) 
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which are characterised by the two mutually orthogonal unit vectors of de- 
formation 

e^ 1 - 1 = e x cos(2na + 9q) — &y sin(2na + 60), (45) 
eff = e x sin(2na + 9q) + e y cos(2na + 60), (46) 

and the amplitudes A x = y / 2/(MN)Q relta , A 2 = \Z2/(MN)Q iml>a . The 
interaction of the quasiparticles with these normal coordinates is determined 
by the coupling functions 

a 

F 2 (k,0) = 4sin(-)( X i + G 2 cos(fc)), (47) 

1— CX OL Oi 

Fi(fe,±a) = ±2V2isin(a)[xicos(-) +x 2 sin(-) + G 2 cos(/c± -)]. 

(48) 



For other frequencies with q > a we have 

^1,2 (Q) = \(^i+ "2 T V /A2 + 4 / 2 ) . ( 49 ) 
and the electron-phonon coupling functions Fj (see (38)) in which 



1 A 

cos(g( g )) = — + — (50) 

siu(%)) = (51) 
/A2 + 4/2 + A)v/A 2 + 4/ 2 

4 Equations in the adiabatic approximation 

Next we consider the ground state of the circle with one quasiparticle. We as- 
sume that the adiabatic approximation is applicable and that we can choose 
the wave-function of the ground state in the form 

|V) = U\4> e ), (52) 

where U is the unitary operator of the coherent molecule displacements, 

\Mt)) = J2Mt)Bl\0) =J2Hk,t)Bl\0), (53) 

n k 
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where the functions (j) n and ip(k) satisfy the normalisation condition: 

]>>ni 2 = i, 5>(*)i 2 = 1 - (54) 

n k 

Here the functions 4> n and ip(k) are connected by the same unitary trans- 
formation from the site representation to the quasimomentum one, (23), as 
the operators B n and Bj,. 

The zero-order adiabatic approximation is equivalent to the semi-classical 
approach when the molecule displacements and their canonically conjugate 
momenta are considered as classical ones. Indeed, by computing the mean 
value of the Hamiltonian over the state (52), which is a functional of (fr n , and 
of the molecule displacements and of their canonically conjugate momenta, 
we transform the problem to a classical one which involves classical variables 
instead of operators. In the site representation the Hamiltonian functional 
then reads as 

n 

+ G 2 {ip*ip i+1 + ip* +l i Pi )U n ) + W, (55) 

where all the variables are c-numbers and W is the mean value of the phonon 
Hamiltonian, H p h, i.e. of the classical Hamiltonian of harmonic vibrations 
of the circle. We can also represent this Hamiltonian in the quasimomentum 
representation 



k j,q 

+ - 7 ^T,^( k ^W(k + q)m(P q , 3 +P- q , j ), (56) 

where the variables tp(k) and {3 q j are connected to (fr n and the molecule 
displacements U{ and Sj and their canonically conjugate momenta p n and q n 
by the same canonical unitary transformations (23), (24), (29), (34) as the 
corresponding operators. 

Then we can write down the Hamiltonian equations of motion for these 
variables. In particular, we get 

" - M ir- 
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We will restrict our attention to the stationary solutions for which (j) n {t) = 
ex.-p(—iEt/h)4> n (or ij)(k,t) = exp(—iEt /h)ijj(k)) and p n = q n = 0. Perform- 
ing the variations of (55) with respect to the variables (fi, Ui and Sj, we 
obtain the static equations 

= (£ - E)tpi - J {tp i+1 + ipi-i) + xWi {Ui + Ui-i) +X2fiAi 

+ G 2 {vi+iUi + ipi-xUi-x), (58) 
= k(cos 2 (^)(2^ - u i+ i - m-i) + cos(-)sin(^)(sj_i - s i+1 )) 

+ K 6 (sin 2 (^)(2ni - u i+2 - Ui- 2 ) + cos(^)sin(^)(2si_i - 2s i+ i + s i+2 - Sj_ 2 )) 

+ (Xicos(-) + X2sin(-)) (|^i-i| 2 - |<^+i| 2 ) 

+ G 2 cos(-) (ift-ifi + <Pi<Pi-i ~ <Pi+i<Pi ~ <Pi<Pi+i), (59) 

= Ksin(-)(sin(-)(2si + s i+1 + Sj_i) + cos(-) (u i+1 - iti-i)) 

+ K b [cos 2 (^)(6sj - 4s i+ i - Asi-i + s i+2 + Sj_ 2 ) 
+ cos(-) sin(-)(2«i+i - 2m-i + m- 2 - u i+2 )\ 

+ 2 (xisin(|) +X2Cos(|))|^| 2 + (xisin(|) - X 2 cos(^))(|i^+i| 2 + |^_i| 2 ) 

+ G 2 sin(— )(^*_ 1 < / o i + p-Vi-i + Vi+iVi + ¥>*¥>*+i)- (60) 

Note that the equations (59-60) are invariant under the following transfor- 
mation: 

Hi > Ui, 
Si ► Si, 

Xn -> -Xn (n = 1,2), 

G 2 -► -G 2 . (61) 

In the quasimomentum representation this system of equations can be 
rewritten in the following form 



(£(k) - E) ip(k) + ]T F j( k ~ 9, q)Qj(q)W " <?) = 0, (62) 

1,3 

tf{q)QM = E F i iWikMk + q), (63) 
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where 

QM = J^(f> q j + (r- q j)- (64) 

From Eq. (63) we obtain the explicit expressions for Qj{q), namely: 

w=-^x;2gV(w(*+ 4 ). (65) 

Substituting these expressions into (62), we obtain the following nonlin- 
ear integral equation for tp(k) 

(S(k)-E)^(k) -± ]T GfakuqWihMh+qMk-q) 
g,fci(|g|>a) 

fcl 

- ^ E ^(*.*i.?)f(*i)#i+?)#-?)=0. (66) 

fci,g=±a 



where 



G(k,h,q) = £ , (67) 



F 2 ( fc,0)F 2 *(fci ,0) 
Aft2|(0) 



MKh) = ^'T/rtZ' > (6 8 ) 



?a(k,h,a) = j^T^ • (69) 

Taking into account expressions (32) and (38), and then Eqs.(28) and (39), 
we obtain 

G(k,k 1 ,q) = ^+G 1 {k,k 1 ,q), (70) 

G 1 (k,k 1 ,q) = -{ x lcos 2 (h + X iG 2 [cos(k - f ) + cos(fci + §)] (71) 
k z z z 

+G\ cos(/c - |) cos(A;i + |)}. 
Substituting expressions (43), (44) and (48) into (68) and (69), we get 
F (k,h) = Gi(Mi,0), 

4 a 

F a {k,ki,a) = Gi(k.ki,a) + -sm(-)f(k,ki,a), (72) 

K 4 
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where f(k, k±,) is a function of k and a, the explicit expression of which will 
be not essential for us. 



5 Fully delocalised states 

Equations (59-60) or (62-63) always admit a simple solution for all values of 
the parameters that correspond to a fully delocalised quasiparticle state, i.e. 
for which the quasiparticle is evenly distributed along all the chain: ipi = 
l/(y/N) exp(ifcon). In the /c-representation this solution is fully localised 
ip(k) = 5k,k - From (63) we see that in this case Qj(q) ~ <5g,o and, therefore, 
only the normal coordinate (22(0) is nonzero. For j > the lowest energy 
state corresponds to the value fco = 0. Defining 

A = E - £ + 2 J, (73) 

this fully delocalised solution for the ground state is given by 

i>{k) = S kfi 

O m 1 ^(0,0) = { [m X i + g 2 

^ 2[ > y/MN fil(0) " V^srn(f) {<V 

in the /c-representation, and 



ATl/2' 

Ui = 0, 



1 Xi + G 2 
Nk sin ^ 



(75) 



in the site representation with 



A = - 4te + G2)2 =AJ. (76) 

According to (73), A = A J determines the eigenvalue of Eq. (59) (or (62)) 

E = £ - 2J + A = 8(0) + A, (77) 

where £ (0) = £q — 2J is the lowest energy level (27) of a quasiparticle in 
the circle. 
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The total energy of the system is Etot = E+W. For the fully delocalised 
solution we have W = (1/2) ^ q ^]{q)\Qj{q)\ 2 = ( 1/2) 0| (0) | Q 2 (0) | 2 . So, 
the total energy of the delocalised state is 

Ekot = £o ~ 2 J - 2(Xl A ^ 2)2 = f (0) - eJ, (78) 

where e characterises the decrease of energy due to the electron-phonon 
interaction. Thus, we see that the localised state of a quasiparticle in a 
circle is characterised by the following values of A and e 

\ 2 9l 9l / vn x 

Mel = > e del = Jj, (ly) 

with 

9i = 2to + «! (80) 

In this state all the s n are the same and all the u n vanish; i.e. the 
circle radius changes due to the alteration of the distance between sites 
U n = —jj^ixi + G2). Note that, in the limit of large N, sina pa 2n/N and 
sin 2f « 7r/AT. Thus, Sj = — (xi + G^VC^) which is independent of q. 



6 Numerical Solutions 

The discrete system of equations (59-60) can be solved numerically. As 
mentioned above, the system always possesses a fully delocalised solution in 
which the quasiparticle is evenly distributed between the chain sites. But 
for most values of the parameters, this is not the solution with the lowest 
energy. In this section, we describe the numerical solutions (we have fixed 
the following values J = 1, k = 1 and kj, = 0.1). 

We have found several regimes of quasiparticle localisation. The quasi- 
particle can be strongly localised on one or more sites (Fig. la), delocalised 
over most of the chain like (Fig lc). Moreover, the maximum of \(fi\ 2 can 
be located on one site (Fig la) or two sites (Fig lb). Like all the solutions 
that we have found, these solutions have a reflection symmetry, but in the 
first case the axis of symmetry passes through a circle site while in the other 
one, it passes in between two sites. 

We would like to point out here that the actual position of the soliton 
on the chain is completely arbitrary and on each of the solution presented 
here it was picked up more or less randomly by the relaxation method. 



17 



Moreover, one should note that given our choice of signs (for its definition), 
A corresponds to the binding energy and so the ground state is the one for 
which A is the largest. 




c d 



Figure 1: Numerical solutions of Eqs. (59-60) for G2 = 1, X2 = 0.2 and 
a) xi = — 2 b) xi = —1-5 c) xi = —0.5 d) xi = 2. Symbols + represent 
the undeformed chain. Symbols in the form of small circles represent the 
deformed chain, while the continuous line represents \(fi\ 2 plotted on the 
undeformed chain. In d), we also show where each site has been moved to 
by the chain deformation. 

One also notices that the chain can be deformed in different ways. For 
example, in Figs, la and lb, the chain is pinched into the shape of a bean 
while in Fig. Id the chain is compressed with a kink appearing where the 
soliton is located. The deformation of the chain in all the figures presented 
here has been scaled in exactly the same way to emphasise the deforma- 
tion of the chain and show the relative magnitude of the deformation when 
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the model parameters are varied. The actual amplitude of deformation de- 
pends on the values of the parameters and is much smaller than the distance 
between the chain sites. 

Notice also that by virtue of the symmetry (61), the solution in Fig. l.a 
can be transformed into a solution for G2 = — 1, X2 = —0.2 and Xi = 2 
where, as the sign of u and s changes, the chain is pinched on the sides so 
that the lattice sites where the soliton is located move outside the circle. By 
the same symmetry, the solution shown in Fig l.d is transformed into a ring 
that is stretched out. 

For large values of xi and X2, the solution is usually strongly localised 
while in the region xi = —G2, the soliton is fully delocalised. This is shown 
in Figs. 2-4 where we represent the value of the maximum of |^| 2 as a 
function of xi and X2 for different values of Gi- 

When G2 = 0, Fig 2, one can see that if X2 > 0.55, the soliton is always 
strongly localised, while for smaller values, the soliton is strongly localised 
when xi << 0. As xi increases, the quasiparticle becomes more and more 
localised, then it becomes fully delocalised around xi = —0.5. When xi 
is larger than 0.5, the quasiparticle becomes localised again and it is fully 
localised as xi g° es to 00. 

At G2 = 0.5, Fig. 3, we see a similar pattern as for G2 = except that 
the region where the soliton is delocalised is centred around xi = —G2- For 
small value of X2, there is also a region where the quasipartivle is delocalised 
(—1.5 < xi < —1 when X2 ~ 0). The solutions in that region are delocalised 
but have the same symmetry as the fully localised solitons, i.e. the axis of 
symmetry passes through a chain site. In the positive range of xi there is 
also a region where the quasiparticle is delocalised (0 < %i < 1 when X2 ~ 0) 
but in this case the solution has a different symmetry, i.e. \ifi\ 2 takes its 
maximum value on two sites and the axis of symmetry passes between two 
chain sites. 

When G2 = 1, Fig 4, we can distinguish 4 different regions. On the 
outside, the black region corresponds to fully localised solutions (Fig. la). 
On the left, there is a sharp grey region (— 1.8 < xi < 1 at X2 ~ 0) where 
the quasiparticle is mostly localised on two neighbouring sites (Fig. lb). In 
the white region, the quasiparticle is fully delocalised (Fig. lc) while in the 
grey half circle, on the right hand side, the quasiparticle is delocalised with 
the axis of symmetry passing between two chain sites . Outside this circle 
the quasiparticle is localised on one site, hence, its symmetry is different 
(Fig. Id). 

When G2 = 1.5, the structure of the solutions is similar to that of 
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G2 = 1, but it now exhibits a region where more than one solution exists 
(other than the fully delocalised solution that always exists). In Fig. 5 we 
plot the value of \(fi\ 2 as well as A as a function of xi an d one notices an 
even richer structure. Notice that A corresponds to the eigenenergy with 
opposite sign, so a large value of A corresponds to a lower energy state. 

In Fig. 6 we present several solutions for the same value of G2 = 1.5. 
When xi = — 2.7, Fig.6a, \(fi\ 2 = 0.803 and the quasiparticle is strongly 
localised on one site. When xi = —2.65, Fig. 6b, \<fi\ 2 = 0.50 and the 
quasiparticle is strongly localised on two neighbouring sites. When xi = 
—0.95, Fig. 6c and Fig. 6d, there are two solutions, \(fi\ 2 = 0.49 and 
\(fi\ 2 = 0.13. The first one is localised on two sites, while the second one 
is delocalised and has a different symmetry. Fig. 6b shows that the most 
localised of these solutions has the lowest energy. When xi = 0.05, |y?i| 2 = 
0.43, Fig. 6e, ie just before the little spike in Fig. 5a, the quasiparticle is 
slightly delocalised and has a symmetry axis passing through a chain site. 
When xi = 0-1) Wi\ 2 = 0.38, Fig. 6f, ie just after the little spike in Fig. 
5a, the quasiparticle is still slightly delocalised but it possesses a different 
symmetry, ie the axis of symmetry passes between the 2 chain sites. When 
Xi = 2.7, \ipi\ 2 = 0.48, Fig. 6g, the quasiparticle is strongly localised on two 
sites. When xi = 2.75, \ipi\ 2 = 0.714, Fig. 6h, the quasiparticle is localised 
on one site. 

In Fig. 7 we present the value of | A| 2 as a function of G2 and X2 for Xi = 
while Fig. 8 presents some typical solutions for that range of parameters. 
We notice that when X2 and G2 are small, the quasiparticle is always fully 
delocalised. When X2 is larger than 0.5, there is a triangular region where 
the solution is strongly localised on one lattice site. We note also that there 
are 3 different zones with the same shade of grey which we have marked as 
III, II and 3 in Fig. 7 and which correspond to solutions spread mostly on 
3 sites (Fig. 8a and 8b), 2 sites (Fig. 8c and 8d) and 3 sites (Fig. 8e and 
8f), respectively. 

The most surprising solutions, shown in Fig. 8g and 8h, are non symmet- 
ric. They occur when X2 > 0.5 and when G2 is large. When G2 decreases 
for a fixed value of X2 , the asymmetry of the solution decreases so that the 
solution becomes fully symmetric and the quasiparticle is localised on two 
sites. 

Our numerical solutions show that there are various types of solutions. 
In particular, they show that if %i, X2 and G2 are not too small, the electron- 
phonon interaction is responsible for the localisation of the quasiparticle. In 
particular, a strong interaction between the quasiparticle and the phonon 
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Figure 5: Solutions of Eqs. (59-60) at G 2 = 1.5 and \2 = 0.2. The plots of: 
a) b) A. 
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Figure 6: Numerical solutions of Eqs. (59-60) for G2 = 1.5, X2 = 0.2 and 
a) xi = -2.7 b) xi = -2.65 c) X i = -0-95 (A = 3.38) d) X i = -0-95 
(A = 0.09) e) xi = 0.05 f) xi = 0.1 g) Xi = 2.7 h) X i = 2.75 . The + 
represent the undeformed chain. Symbols are explained in caption of Fig. 
1. In (g) and (h), we also show where each site has been moved to by the 
chain deformation. 
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j.1 M kb.0.1X1.0 



X2 




I 0.998783 
I 0.903905 
I 0.809026 
I 0.714148 
I 0.61927 
0.524392 
0.429513 
0.334635 
0.239757 
0.144878 
0.0500001 



Figure 7: Density for the solutions of Eqs. (59-60), |</?j| 2 , for xi = 0. 



fields {G2 large) introduces some sharp transition between the different types 
of solutions and also leads to an overlap between them. 

It is also interesting to note that overall, the solution is either localised 
on a few lattice sites or spread nearly uniformly across the entire chain, 
but the range of parameters for which a broad soliton exists is very narrow. 
This is shown on Fig. 2 -4 by the sharp transitions between the grey and 
the white zones. 



7 The Continuum Limit. 

Though the system of equations always admit localised solutions, these so- 
lutions are not always the states of the lowest energy. At strong enough 
electron-phonon coupling the self-trapping of a quasiparticle can take place 
with the formation of polaron-type states. Depending on the strength of the 
coupling, the localisation can be (i) strong with the formation of a small po- 
laron state with the localisation of a quasiparticle mainly on one lattice site, 
or (ii) comparatively weak with the formation of a large polaron state with 
the distribution of self-trapped quasiparticle among several lattice sites. 
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Figure 8: Numerical solutions of Eqs. (59-60) for \i = > a ) X2 = 0.1, 
G 2 = -3 b) X 2 = 0.1, G 2 = 3 c) X 2 = 0.5, G 2 = -2 d) X 2 = 0.5, G 2 = 2 
e) X 2 = 0.5, G 2 = -0.7 f) X 2 = 0.5, G 2 = 0.7 g) X 2 = 0.5, G 2 = -7 h) 
X2 = 0.5, G 2 = 7 Symbols are explained in caption of Fig. 1. In (g) and 
(h), we also show to where each site has been moved to by the deformation 
of the chain. 
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We consider this latter case here. In linear chains the solutions are well 
described in the continuum approximation which is valid for large enough 
systems i.e. N 3> 1. For a circular chain this means also that a is very 
small. For such circles, with large enough N, we can neglect the term 
involving sin(^) in (72). In this case Eq. (66) can be rewritten in the form 

(£{k) E)^{k) -^Giftti.^^ + ^i-?) 

q,ki 

^ E -#i)#i+9)#-?) = 0. (81) 

g,fci(|g|>a) 

In the last term on the r.h.s. of Eq. (81) we can use the approximation 
^2q,ki,\q\>a ■ ■ ■ = ^q,ki ■ ■ ■■ Of course, the value of the constant b 2 <C 1 
depends on a as well as on the degree of localisation of the solution. For 
fully delocalised solutions b 2 = but b 2 rapidly tends to unity as the solution 
becomes more localised. In this way Eq. (81) becomes 

(£(k)-E)^(k) -l-Y^GikM.qWikMki+qmk-q) = 0, (82) 

q,ki 

where 

2 

G(k,h,q) =b 2 — + G l (k,k 1 ,q). (83) 

K b 

A weak localisation in the site representation corresponds to a strong 
localisation in the ^-representation, i.e. ip(k) (k = av) is essentially nonzero 
only for several values of u in the vicinity of v = 0. In the long- wave 
approximation when k = av <C 1 (which assumes N S> 1) we can write 
down the following expansions for the functions in (66): 

S(k) = £ - 2J + Jk 2 + . . . , 

G(k, h,q) = b 2 ^+4 iXl+G2)2 +.... (84) 

Kb K 

To solve Eq.(66) we introduce the function 




of the continuum variable x. Note that ip(x) is a periodic function, <p(x + 
N) = ip(x) and, at x = n, this is a unitary transformation of ^^{k) to the 
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site representation. Moreover, if <p(x) is known, we can find ip(k). 



i 2 



^(fc) = e- ikx <p(x)dx. (86) 



Using the approximation (84), one can transform (66) into a nonlinear 
differential equation for ip(x): 

jd ~^T + G \^)\ 2 v(x) + M*) = 0, (87) 

which is the stationary nonlinear Schrodinger equation (SNLSE). Here A is 
defined in (73) and 

G = G(0, 0, 0) = 4 (X1 + + (88) 
Next we rewrite (87) in the dimensionless form after defining 
A = -n J, and G = 2 J5 where g = gi + g b . g b = ——, (89) 
where gi is given in (80). We get 

^fi + 2 5 | ¥ >(x)| - = 0. (90) 

We look for a solution of (90) which satisfies the periodic boundary- 
condition 

ip(x + L) = <p(x), (91) 
and the normalisation condition 

[ L \<p(x)\ 2 dx = l, (92) 
Jo 

where L is the length of the ring (in the adimensional units used here, 
L = N). We see that ip(x) is a real function and so we can integrate Eq.(90) 
and obtain 

<dip(x)\2 
dx 

where C is the integration constant. So we have to invert 



ra=,vw-^)+c (93) 



= _ r (x) ^ 04) 

Mo) vW- W 4 + C 
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Using the substitution 



V 2 = -4~T), (95) 



<T 3 

we reduce the problem (94) to a normal Weierstrass form 

f^(§) dr 



['w dr , . 

/ / 5 = ( 96 ) 

•ArM v 4r J - qot — Q?. 



'HO V 4 T 3 - 52T - 53 

which is solved by the Weierstrass elliptic function p(x + c) [20, 21] (c 
being some constant). Here, and in what follows, we shall use the standard 
notations for the theory of elliptic functions. In our case, the constants g 2 
and 53 in (96) are determined by 

92 = 4(i M 4 + 0C), 

53 = -\f{\i+ + gC). (97) 

The elliptic Weierstrass function p(z; 52, 53) = p(z\w, a/) is a two-parametric 
doubly periodic function in the complex plane z = x + iy. The two param- 
eters are either the constants 52 and 53 or, equivalently, the two periods, 
which are denoted as 2ui and 2a/, of the function p(z + 2mw + 2na/) = p(z) 
with m and n being integers. In the theory of elliptic functions [20, 21] the 
following symmetric notations are generally used 

uj\ = lo, UJ2 = —uj — a/, CJ3 = a/ (98) 

so that p{uij) = ej (j = 1, 2, 3) where e,- are the roots of the cubic equation: 

4r 3 - g 2 T - g 3 = 0. (99) 

By definition, \ip\ 2 is a real and bounded function of the real coordinate 
x. Taking into account the analytical properties of the Weierstrass function, 
we conclude that 99 defined by (96) will be physically meaningful only when 
the discriminant g% — 27^3 of the cubic equation (99) is positive and all three 
roots ej are real and different. In the theory of elliptic functions [21] the 
real roots ej of the cubic equation (99) are distributed as follows: 

ei > e 2 > e 3 and e_? > 0, e 3 < 0. (100) 
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In our case g\ — TIg\ = IQg 2 C 2 (fj, +4gC) and the roots of the cubic equation 
(99) are 

i(-i/x 2 ±vW4 5 C), (101) 

i.e. for the integration constant C we have the constraint AgC > —/J, 2 . In 
this case the period 2u> is real, 2a/ is purely imaginary, and the solution of 
(96) reads as p(x + a/). Therefore the solution of Eq.(93) (or (94)) is given 
by 

¥> 2 (*) = -(^-p(*W)). (102) 
9 

This solution must satisfy the boundary (91) and normalisation (92) condi- 
tions which lead to the relations 

2mw = TV, 
-H 2 N + 2mri = g. (103) 

Here m = 1,2,... is an arbitrary number and rj = ((uj) is the constant of 
the elliptic functions theory (((z) is the Weierstrass ("-function determined 
by the relation p(z) = —Q'{z)). The conditions (103) allow us to determine 
the integration constant C and the eigenvalue fj, 2 , and then to calculate the 
total energy of the system which includes the energy of the deformation. 
Using solution (102) we get the following expression for the total energy: 

E tot = £ (0) - eJ, e = ^ 2 - l -CN. (104) 
To solve the system of equations (103), let us introduce the parameter 

k 2 = ^Lll < i (105) 

ei - ea 

which, from now on, will be considered as an independent parameter instead 
of the integration constant C. Here e/ = p(tOj) are the roots the cubic 
equation (99) and k is the modulus of the elliptic Jacobi functions [21]. 
Using the formulas 

K(k) = ( ei - e s )iu, E(fc) = + (106) 

(ei - e 3 )z 
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where K(fc) and E(fc) are complete elliptic integrals of the first and second 
kind, respectively, we can rewrite the conditions (103) in the form 

2mK(k) = N( ei - e s )i , 
2m(e 1 -e s )*E(k)-(e 1 -^ f i S )N = g. (107) 

Inserting the first expression into the second one, we can rewrite (107) as 

2mK(k) = N( ei - e 5 )i, 
(e, - e s )E(k) - (e, - |// 2 )K(fc) = ^-( ej - e 5 )i (108) 

The solution (102) can be rewritten in the form 

^) = !((£-„)- ("- ")("-") ) . (109) 
g\ 3 p(x)-e 3 > 

Depending on the sign of the integration constant C, the roots (101) will 
be enumerated in two different ways according to the rule (100): 



1 2 



62 = r 



-1(^4 + 4^+1 2 } (no) 



at C > and 



1 2 

^ = 3M, 

e 5 = -^(^W + / 2 ) (HI) 

at C < 0. Two different solutions correspond to these two cases of root 
distribution. 

1) C > 0. In this case, taking into account the root distribution (110), 
we obtain from (105) the following expression for the constant C in terms 
of k: 

C ~ {2k* -1)* J' { } 
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The second condition in (108) gives us the expression for \i 

g ^/2k 2 - 1 

/i ~ 2^B{k)-(l -k 2 )K(k) (U3) 

and the first one gives us the equation which determines the constant k, the 
values of which are restricted by the condition that k > jj , namely: 

K(i)(E(fc)-(I-i')K(t)) = i (114) 

Using the expression for the elliptic Jacobi functions in terms of the 
Weierstrass function, 

cn 2 (U; k) = P( x ) ~ e i u = ( ei - e3 y/2 x (H5) 
p{x) - e 3 



we find that the solution (109) is given by 

"2mKx 



N 



,k 



(116) 



From the periodic properties of the function cn(u,k), cn{u + 4K,&) = 
cn(u, k), it follows that the periodic condition (91) will be satisfied if m = 
2n, n = 1,2, .... The total energy (104) of the nanocircle in these states is 
then 

E tot (n) - £(0) - — ^ _____ . (117) 

2) C < 0. In this case the roots are distributed as (111) and from (105) 
it follows that 

1 - k 2 ii 4 

c -jh^rT < 118 » 

The second condition in (108) gives us the expression for fi 



g V2 - k 2 , s 

and the first one gives the equation for determining the constant k 

E(t)K(t) " Wf- (120) 
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Taking into account the relation 

dn 2 (u,k)= pi / X) ~ e2 , (121) 
P(x) - e 3 

we find that in this case a normalised periodic solution can be expressed in 
terms of the elliptic Jacobi function as follows: 

~2mK(k)x 



N 



(122) 



Using the periodic properties of the function dn(u,k), dn(u + 2K, k) = 
dn(u, k), it follows that the number m can take the values m = 1, 2, .... The 
total energy (104) of the nanocircle in these states is then 

p ( m ) — £(q) J 9 2 E(AQ(g-fc«)+K(fc)(i-fc«) 

E tot (m) - £(0) - ^-j . (123) 

The analysis of Eqs. (117) and (123) shows that the energy of states 
(122) is always lower than the energy of states (116). Because we are inter- 
ested in the solutions which correspond to the lowest energy we can ignore 
the solution (116) and consider only the solution (122). Moreover, we are 
interested only in the ground state of the circle which corresponds to the 
value m = 1 in (122), (123) and (120). All other states are excited states of 
the circle. 

From (119) and (123) we obtain the following expressions for the eigen- 
value 

-"T^p (124) 



and for the total energy 

g 2 E(k)(2-k 2 )+K{k)(l-k 2 ) 



(125) 



12 E 3 (k) 

of the quasiparticle ground state in the chain. 

Let us consider next the case when the modulus of the elliptic function 
is small, k 2 <C 1 , which corresponds to the case of small nonlinearities and 
not too long chains. In this case, using the expansions of complete elliptic 
integrals in series of k, 

KW = f(l + l*' + ^-...), 
E(*) = !(l-1*«-A*4_...) (126) 
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we have from (120), at m = 1 that 



k*=tV20-l. (127) 

So we can conclude that the solution exists only when g exceeds the critical 
value of the nonlinearity constant 

9cr(N) = ^ (128) 

which depends on the length of the ring. 

In this case the solution (122) can be written as 



t(i-v / ¥^ 



V' 2 ^(i-v^) 

(129) 

For the eigenvalue and the total energy at 2 <C 1 expressions (124) and 
(125) give 



A^-^l + i^), ^^(l-JLjfc*), (130) 

Remember, that this solution exists only when g exceeds the critical value 
g cr (128). Below g cr the fully delocalised state is realised on the circle. 
At g — > g cr the solution (129) formally tends to a delocalised solution but 
comparing (79) with (130) shows that even at g = g cr the total energy of 
the state (122) is lower than energy of the fully delocalised state. Only in 
the case = does the energy of the fully delocalised solution coincide 
with the energy of the periodic solution at g = g cr . At g > g cr the energy of 
the periodic solution is always lower than the energy of the fully delocalised 
state. 

Let us consider now the opposite case, i.e. when the modulus of the 
elliptic function is close to 1, k 2 « 1 , kf = 1 — k 2 <C 1 . This case, as we 
will see below, corresponds to long chains and not too small nonlinearity 
parameters, such that the condition giV » 1 is satisfied. Then, using the 
expansion of the elliptic functions and integrals in power series of kf , 

K(k) = In ^- + ^ fin j- - 1 )kf + . . . , 



ki 4 \ hi 

l( ln ^ 2 
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we get from (120): 

k, = 4e~ 9N / 4 , k 2 = 1- 16e~ 9N / 2 . (132) 
Substituting this into (122), we get 

^ * 2[l + 2( 5 iV^2) e -^/2] [1+ 



_l_ e 5^/4 — . ^ 9_ Yi — ) tanh — 
V 2 2 2 / 2 



cosh" 1 ^, (133) 
For the eigenenergy and the total energy we get from (124) and (125) 

A d „«-^(l-8(2 5 iV-3)e-^/ 2 ), e dn « g (l + 24 e -^/ 2 ) . (134) 

We see that at large enough (7 (giV 3> 1) the solution (116) describes 
a well localised state and that increasing g leads to the increase of the 
localisation. But it is necessary to keep in mind that the applicability of the 
continuum approximation is restricted by the condition that the localisation 
is not too strong which itself imposes the constraint \i < it/ 2. 

Using the obtained solution it is possible to calculate the configuration 
of the circle in the state (116). Taking into account (24), (29), (65), and 
the explicit expressions for the corresponding functions (32), (38), (39), (28) 
and (30) we can write 

l^e^t/^), s n = l^e^), (135) 



— - N q 



with 



S in(f ) cos(f ) + cos(f ) sin 2 (§ )[£ cos(f ) + ^ cos(fe + §)] 
sin(f)[sin 2 (f)cos2(f)-cos2(f)sin 2 (f)] 



4k 6 



U ^ q) % ^ -'-^m 2 (§)cos 2 (f)-cos 2 (§)sin 2 (§ 

xip*(k)ij(k + q), (136) 



^ cos(f ) - sin(^) cos(f) [2£ cos(f ) + cos(A: + §)] 

(137) 

Here VK^O is given by (86) in which ip(x), in turn, is given by (122). 

In Fig. 9 we present the curve of g = g cr = 7r 2 /20 for N = 20, (89), 
derived in the continuum approximation, where G is given by (84) with 
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-2 -1 1 2 -4-3-2-101234 

XI G2 

a b 

Figure 9: Regions of solutions: delocalised (0), strongly localised (I) and 
localised on two adjacent sites (II) for G2 = 0.5 (a) and xi = 0-5 (b) 
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b = 0.5. The curves on Fig. 9 are the two circles of radius 0.5, but which 
look like vertical ellipses given the different scales of x and y in our plots, 
enclosing the region marked as '0'. All the solutions inside these circles 
(marked as in Fig 9) are delocalised. Fig. 9. a corresponds to the case 
G2 = 0.5 and so can be compared to Fig. 3 obtained numerically, while Fig 
9.b corresponds to the case xi = - to be compared with Fig. 7 



8 The Limit of Strong Localisation 

As we have seen from the discussion above, the increase of the coupling 
constant reduces the region of localisation of the quasiparticle. Thus, at a 
large enough value of the coupling the quasiparticle can be localised mostly 
on one, two, to three lattice sites. In such a case the conditions of the 
applicability of the continuum approximation are not satisfied. To consider 
the limit of strong localization, we use Eqs. (71) -(72) to rewrite Eq.(66) as 

(£(k) - E + 1^! ) m + _L ^[p( a )^(k -a) + p*{a)^{k + a)] 

iV Kfy TV Kb 

^E G ( fc ' fc i'#*( fe iM fc i + ?)#-?) = °- ( 138 ) 

q,ki 

Here, assuming a large enough circle with N ^> 1, we have neglected the term 
~ sin(a/4) in (72). Moreover, we have taken into account the normalization 
condition (54) and introduced the notation 

p(q) = J2r(k)iP(k + q). (139) 



The function ip(k) = (1/yN) J2 n ex P ( — ikn)(f) n is periodic in k, ip(k + 
2tt) =tp(k), and so it is representable as a Fourier series 

oo oo 

ip(k) = -j={cq + ^2 c m cos(mk) + i ^ s m sin(mfc)). (140) 

m=l m=l 

For p(q) we then have 

^ oo oo 

P(q) = c o + 7; ( c « + s ™) cos ( m Q) + i c m s m sin(mq) = p r (q) + ipi(q). 



2 

m=l m=l 



(141) 
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The coeficients of the expansion satisfy the normalisation condition: 

]T \m\ 2 = P(0) = c§ + - £ + 4) = 1. (142) 

k m=l 

Next we put (140) into (138) and get the Fourier series 

{F(0) + ^F(m)cos(mA;) + ^ ip(m) sin(mfc)} = 0, (143) 

m m 

thus showing that F(m) = tp(n) = 0. These conditions give us the equations 
for c m and s n . We have 

A c - J+(0)ci - T_(0)si = 0, (144) 
Aisi - \j+{2)s 2 - N lCl - G c (0)c - \t-{2)c 2 = 0, (145) 
Aici - J c (0)c - X -c 2 - N lSl - ^T_(2)s 2 = 0, (146) 

A m C m - ^J_(m- l)Cm-i - ^J+(m + l)c m+ i - 7V" m Sm 

- ir+(m- l)a ro _i - ^T_(m + l)s m+1 = 0, (147) 

A m s m - -J_(m -l)s m _i - ^J + (m + l)s m+ i - iVmCm 

- iT+(m-l)c m _i - ^(m+l)c m+1 = 0. (148) 

The functions A fc , J±(i), T±(j), N m , J c (0) and G c (0) are functions of s m 
and c n and so we see that we have nonlinear equations for Sk and Cj. 

Note that when we have a fully delocalised system all are the same; the 
same is true of Sk- As the coupling constant increases, s n become different 
from each other; in fact, as n increases, they become exponentially small (as 
a function of n). The same is true also for c n . Thus, to get an approximate 
solution, we can take into account only a finite number of Sj and Cj. 

Our system of equations (146) always allows us to put s n = and c„ / 0. 
If we now assume that c 2 <C c\ < cq, we can consider the equations involving 
only Co and c±. In this case we have 

{A — (— H )cg - — c{ c ci}co 

ft ZK TTK 

- {J + ^ C oc 1 + ^(cg + i)} Cl = 0, (149 ) 

K 7TK D 
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{ 2J+ 2 M CoCl + ^ (c 2 + l )}co = , (150) 

where A = £o — E 1 . 

When ci < Co then Co and Ci are given by 

2 = (i + g 3 )(^-gi+g 2 ) fl51) 

1 D + g 3 (2D + 9l - g 2 y V ; 



{l + Sg 3 )(D + g 1 -g 2 ) 



(152) 



where 5l = + 52 = ^ + _i , 53 = 4*±-f . Here D is given by 



2[L> + 53(2,0 + 5l -g 2 )Y 

X? , ^±1 _ . X1 G 2 
«i,J 1 ftJ ' ^ ~~ kJ ^ k.J ' y 3 ~~ 41 ttkJ 

= v / (5i-52) 2 + 8(l+53)(l + 3<73)- (153) 
The quasiparticle wavefunction is now given by 

ipn = co<5 n , no + ^Ci(5„ j7l0+ i + (5 nj7l0 _i), (154) 

and describes the localisation of the quasiparticle on at most 3 sites . 

The energy of this state is given by E = £ — \{g\ + g 2 + D}. Of course, 
our expression is valid when c\ < c$. 

<? 

In Fig. 9 we present the plots of the curve corresponding to -\ = 0.2 

c o 

obtained from (151), (152) and (153). In Fig. 9. a we present the case 
G 2 = 0.5 as a function of xi an d X2 while Fig. 9.b corresponds to the case 
Xi = as a function of X2 and G 2 . Thus, our plots shows the curves which 
separate the regions of strongly localised solutions from the region where the 
quasiparticle is localised on a few adjacent lattice sites. The curves are the 
ones that delimit the region marked as I in the figure and which corresponds 
to the strongly localised solutions. When we compare Fig. 9. a and Fig. 9.b 
with Fig. 3 and Fig. 7, respectively, we note a very good agreement between 
our numerical results and our analytical results. 

Finally, we consider the case where the quasiparticle is restricted to two 
neighbouring sites. In this case we consider equations for cq, c\ and s\ 
assuming that all other Sj and c% are negligeably small. Thus, c 2 = 0(e 2 ) 
and so is s 2 . 
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In this case the equations are somewhat complicated, so here we present 
only their solutions. They are given by: 



d = c + 0(e 2 ), Sl = c - 0(e 2 ). (155) 

The condition of normalisation gives us Co ~ -j=. The quasiparticle wave- 
function is then given by 

V>n = - 0(e 2 )j (5 ntno + S ntno+1 ). (156) 

The energy expression is now given by 

E = So -J ~ \(9i+92) + 2g 3 . (157) 

Our numerical calculations, see Fig. 6c and d, have shown that at the 
same values of the parameters we can have two different solutions, one lo- 
calised essentially on two sites and the other localised over a larger number 
of sites. These two solutions have different energies. 

The necessary condition for the existence of solutions localised on two 
adjacent lattice sites is that g\ — 52 satisfies the relation 

*-»-!7-$ + S>> a <158) 

The physically relevant solution is, of course, given by the state which has 
the lowest energy. The two site solution has the lowest energy when G2 is 
significantly larger than all other parameters. 

In Fig. 9 we present the curve g<i — g\ = corresponding to the limit 
of a solution localised on two adjacent sites. Fig. 9. a presents the case of 
G2 = 0.5 as a function of xi and X2 and the curve is the ellipse, appearing 
as a circle on our scale, near the symbol 'II'. The solutions localised on two 
sites exist just outside the ellipsoidal region marked as iT. Note also that in 
the region contained in the intersection between region and region II, the 
solutions are delocalised. All this confirms the numerical results displayed 
in Fig. 3. 

Fig. 9.b presents the case xi = as a function of \2 and G2. The curves 
are the two straight lines and the solutions localised on two sites exist in 
the region marked as //, just below the two lines. This is also in very good 
agreement with Fig. 7. 
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Notice also that we do not have any theoretical prediction for the tran- 
sition from solutions localised on two sites to those localised on three sites. 

For the case G2 = 0, Eq. (158) has no solution and so there are no 
solutions localised on two neighbouring sites in that case, confirming our 
numerical results. 

We should point out that we have also computed the theoretical predic- 
tions for the different types of solutions corresponding to the case G2 = 
and G2 = 1. When we have compared these predictions to the numerical 
results presented in Fig. 2 and Fig. 4 we have also found a very good 
agreement. 

9 Conclusions 

We have shown that one can successfully describe the electron-phonon in- 
teraction of a circular chain while taking into account the transversal and 
longitudinal displacement of the ring of monomers. To achieve this, we have 
introduced a term in the phonon Hamiltonian which describes the bending 
of the chain. This extra displacement can, in turn, be coupled to the quasi- 
particle and this coupling favours a self-trapping of quasiparticles and the 
formation of localised soliton-like structures. 

The transversal displacement plays an important role in the formation 
of such states. The interaction between the quasiparticle and the phonon 
field leads not only to the localisation of the quasiparticle but also to a 
deformation of the circle. Depending on the coupling between the two fields, 
the circle is pinched at the position of the soliton or on the either side of it. 

We have studied our model by solving its classical equations numerically. 
Our numerical studies have shown a great richness of solutions; hence we 
have also tried to solve our equations analytically. To do this we have had 
to perform several approximations. We have found that, at least in the 
regions of validities of these approximations, our analytical results are in 
good agreement with our numerical results. 

Our results show that, depending of the choice of the parameters, our 
model possesses several classes of solutions. Fully delocalised solutions exist 
for all values of the parameters. They do not correspond to the ground 
state when other solutions exist but they are the only solutions when the 
parameters xi + G2 and \2 are small. On the other hand, strongly lo- 
calised solutions exist when the coupling parameters are large. For these 
solutions, the quasiparticle is localised mostly on one lattice site, but the 
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circle is completely distorted, i.e. the strong localisation applies only to the 
quasiparticle, not the displacement fields. 

We also found some solutions that are localised on two adjacent sites 
and have a different symmetry. 

Our study has shown that the electron-phonon interaction on a circular 
chain leads to the formation of bound states and to the deformation of 
the chain. Our results can be used to the description of various systems 
mentioned in the introduction and as a 'warm-up' stage to the description 
of various 3-dimensional systems such as e.g. nanotubes. 
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